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Abstract 

The neutrino minimal Standard Model (i/MSM) is the minimum extension of the standard model. 
In this model, the Dodelson-Widrow mechanism (DW) produces keV sterile neutrino dark matter 
(DM) and the degenerate GeV heavy Majorana neutrinos lead to leptogenesis. However, the DW 
mechanism has been ruled out by Lyman-a bounds and X-ray constraints. An alternative scenario 
that evades these constraints has been proposed, where the sterile neutrino DM is generated by 
the thermal freeze-in mechanism via a singlet scalar. In this paper, we consider a Higgs sector 
extension of the z^MSM to improve dark matter sectors and leptogenesis scenarios, focusing on the 
thermal freeze-in production mechanism. We discuss various thermal freeze-in scenarios for the 
production of keV-MeV sterile neutrino DM with a singlet scalar, and reinvestigate the Lyman-a 
bounds and the X-ray constraints on the parameter regions. Furthermore, we propose thermal 
freeze-in leptogenesis scenarios in the extended zzMSM. The singlet scalar needs to be TeV scale 
in order to generate the observed DM relic density and baryon number density with the thermal 
freeze-in production mechanism. 
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I. INTRODUCTION 


The standard model of particle physics (SM) has demonstrated great success in high 
energy physics. However, it is not the complete theory because the Standard Model cannot 
treat gravity consistently nor explain several observed phenomena, such as neutrino oscil¬ 
lations, cosmological dark matter, the baryon asymmetry of the universe, the horizon and 
flatness problems, etc. In the dark matter sector, TeV-scale supersymmetry provides weakly 
interacting massive particles (WIMPs) as leading dark matter candidates [1] . However, the 
first run of the LHC experiment excludes a significant region of parameter space for the weak 
super-partners, and the recent results from LUX [2, 3] and XENONIOO [4-6] have severely 
restricted the WIMP cross section. This situation is the same in other beyond standard 
models. Feebly Interacting Massive Particles (FIMPs) [7-10] are not constrained by such 
direct detection experiments due to their much smaller couplings. Therefore, FIMPs are an 
interesting candidate for dark matter. 

A sterile neutrino can be a FIMP. In particular, keV sterile neutrinos are a candidate for 
warm dark matter. The neutrino minimal Standard Model (z/MSM) [11-13] has three right- 
handed neutrinos below the electroweak-scale and the lightest right-handed neutrino may 
have a mass around the keV scale. The lightest right-handed neutrino Ni becomes the keV 
sterile neutrino dark matter and is produced via active-sterile neutrino oscillation, which 
is called the Dodelson-Widrow mechanism [14]. The other heavy right-handed neutrinos 
N 2 and IV 3 lead to leptogenesis via CP-violating oscillations [15, 16], where the heavy right- 
handed neutrinos should satisfy 1 GeV < mjv 23 < 20 GeV. Furthermore, it is possible realise 
Higgs inflation by introducing a non-minimal coupling between the Higgs and gravity [17, 18]. 
Thus, the z/MSM can explain a large number of phenomena with a minimum number of 
parameters. 

However, the z/MSM is severely constrained by recent observations. In particular, there 
are severe constraints on sterile neutrino DM. The Dodelson-Widrow mechanism is known to 
be excluded by Lyman-a bounds and X-ray constraints [19]. To evade these constraints, the 
production of sterile neutrino DM by the thermal freeze-in mechanism has been considered 
in Refs. [20-27] 1 . Ref. [20] considers a scenario in which the inflaton decays into sterile 
1 Ref. [27] gives a comprehensive study of keV sterile neutrino DM via a singlet scalar, but our purpose 
is to estimate the scale in the extended z/MSM to improve the dark matter sectors and the leptogenesis 
scenarios rather than such a generic study of the sterile neutrino DM. 
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neutrino DM. Refs. [21, 22] show that a GeV-scale singlet scalar produces the sterile neutrino 
DM. Refs. [23, 24] consider the non-thermal production via the decay of a singlet scalar. 
Furthermore, the singlet scalar can improve the electroweak vacuum stability or the Higgs 
inflation as well as the dark matter sectors. 

In this paper, we do not discuss the theoretical merits of the singlet scalar in the z/MSM 
such as the electroweak vacuum stability [28, 29], Higgs inflation [30] and scale invari¬ 
ance [31], although we are motivated by these theoretical aspects. Instead, we concentrate 
on estimating the scale of the singlet scalar to improve the dark matter sector and lcpto- 
genesis scenarios. We discuss various thermal freeze-in production scenarios for keV-MeV 
sterile neutrinos in the extended z/MSM with a singlet scalar. In particular, we revisit the 
Lyman-a bounds and the X-ray constraints and show that the singlet scalar cannot be heav¬ 
ier than the TeV scale. We also discuss thermal freeze-in leptogenesis scenarios, which are 
able to produce a larger lepton asymmetry than is produced in thermal leptogenesis due 
to the contribution from the singlet scalar. In the leptogenesis scenarios, the singlet scalar 
needs to be lighter than 1 TeV in order to generate the observed baryon asymmetry. 

This paper is organized as follows. In section II, we review the scalar singlet extension of 
the z/MSM. In section III, we consider two thermal freeze-in scenarios, one utilizing a thermal 
singlet scalar and the other a non-thermal singlet scalar. In section IV, we review the X-ray 
constraints and the lifetime bounds on the sterile neutrino dark matter. In section V, we 
investigate the free streaming horizon and the Lyman-ct constraints in our scenarios. In 
section VI, we discuss thermal freeze-in leptogenesis with the singlet scalar. Section VII is 
devoted to discussion and conclusions. 

II. THE SCALAR SINGLET EXTENSIONS OF THE z/MSM 

In this section, we review the extended z/MSM which contains three right-handed sterile 
neutrinos N a (a = 1,2,3) and one real singlet scalar S 2 . In this model, the vacuum 
expectation value (VEV) of S generates a Majorana mass M a for the right-handed neutrino 
N a . The Lagrangian is given as follows, 

c = C SM + \ (d^S) (d»S) + iW a 0N a - y aa H%N Ra - y SN°N a - V (H, S ) + h.c., (1) 

2 Here we do not consider a complex singlet scalar because in that case a light Nambu-Goldstone boson 
appears with U( 1)l breaking. The presence of such light bosons would make the sterile neutrinos unstable. 
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where i a are the lepton doublets, H is the Higgs doublet, y aa and K a are the Yukawa 
couplings. V (H , S ) is the Higgs potential. After spontaneous symmetry breaking, the Higgs 
doublet and the scalar singlet develop the VEVs (H) = v and (S) = /, respectively, where 
v = 247 GeV and S = s + (S). The right-handed neutrinos acquire the Majorana masses 
Ad a = fta (S). Without loss of generality, we can choose the mass basis where the Majorana 
mass term is diagonal. The Lagrangian is written as follows, 

c = C SM + \ (<V) i^s) + iN^QNR. - y ai H^T a N Ri - ^-N^Nr. - V (H, S ) + h.c.. (2) 


If the Dirac masses are much smaller than the Majorana masses, then, as a result of the 
type I seesaw mechanism, the left-handed neutrino masses can be expressed as follows, 


m u ~ m D M 1 m D ~ , 

k(S) 

mN — Ad~K(S), 6 ~ 


(3) 

(4) 


For the scalar potential V (H, S), we impose the softly broken discrete symmetry Z 2 , where 
the scalar singlet is Z 2 -odd (S —> —S) and all the other fields are Z 2 -even. We can then 
construct the following scalar potential with even powers, 

V (H, S ) = -h 2 h H*H - X -\x 2 s S 2 + A h {H ] H) 2 + ^A 5 ^ 4 + 2A (H*H) S 2 + uS, (5) 


where c uS is a soft Z 2 breaking term. The spontaneous breaking of the discrete symmetries 
Zat could produce domain walls [32], The soft Z 2 breaking term coS makes the vacua of the 
singlet scalar degenerate so that the domain wall problem is evaded [32-34] . The minima of 
the scalar potential are given by the following equations, 


L 2 H = X H v 2 + 2Xf 2 , 
l = X s f 2 + 2Xv 2 + ‘f . 


( 6 ) 


The mass eigenstates of the 
corresponds to the SM higgs 


Higgs and singlet scalar are h and s, where h approximately 
boson. The physical masses of h and s are given by, 

(2A fvf 


m h ~ X H v - A 


m 2 s - A s f 2 + 


As/ 2 — X H v 2 ’ 
(2A fvf 
XsP-Xhv 2 ' 


(7) 

( 8 ) 
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The Higgs portal coupling A induces doublet-singlet mixing. In this paper, we consider a 
TeV-scale singlet scalar which decays into sterile neutrino DM and heavy Majorana neutri¬ 
nos. Therefore, there is essentially no constraint on the coupling A. However, the size of 
the coupling A still affects the thermal history of s. The references [22, 35] show that s is 
out of thermal equilibrium if the Higgs portal coupling satisfies A -C 10~ 6 . In this paper, we 
assume that the singlet scalar is out of thermal equilibrium for A <C 10 -6 and concentrate 
on the thermal freeze-in production mechanism 3 . 

III. STERILE NEUTRINO DARK MATTER FROM THE THERMAL FREEZE- 
IN PRODUCTION MECHANISM 

In this section, we consider various scenarios for the production of sterile neutrino DM in 
the extended z/MSM with a singlet scalar. Sterile neutrinos could be produced by thermal 
freeze-out, thermal freeze-in or non-thermal decay. These scenarios also depend on whether 
the singlet scalar is generated by freeze-out or freeze-in. In addition, the Dodelson-Widrow 
mechanism can produce sterile neutrinos via active-sterile oscillations. 

It is possible to constrain these scenarios using the mass relation of the seesaw mechanism. 
For simplicity, we assume that the lightest right-handed neutrino is sterile neutrino DM, 
with a mass of about 10 keV. We will see later that a sterile neutrino with mass above 1 
MeV is not favored by X-ray constraints and lifetime bounds. If the Yukawa coupling of the 
singlet scalar and the right-handed neutrino is K\ ~ 10~ 8 and the vacuum expectation value 
of the singlet scalar is ( S) ~ 1 TeV, then the following relations can be derived from the 
seesaw mechanism, 


m D M 1 m D ~ ^ — V 2 (lO 18 eV) , 

M {S) 

(9) 

rriN 1 — Mi ( S ) ~ 10 keV. 

(10) 


The Yukawa couplings y, K\ are very small y ~ 10~ 10 and K\ ~ 10 -8 . If the reheating 
temperature Tre satisfies Tre 1$ 10 16 GeV, the sterile neutrino DM does not come into 
thermal equilibrium for K\ <C 10” 6 [36, 37]. Therefore, we may regard the sterile neutrino 
DM as non-thermal particles in the early universe. In such a case, we find only two realistic 
3 If the singlet scalar is not directly produced from inflatons and the reheating temperature is low enough, 
the singlet scalar can be out of thermal equilibrium even if A > 1CT 6 . 
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dark matter scenarios to realize keV-MeV-scale sterile neutrinos. We will now proceed to 
describe these two scenarios. 

A. The singlet scalar is in thermal equilibrium 

If the Higgs portal coupling is relatively large A > 10“ 6 , the singlet scalar s enters into 
thermal equilibrium and the sterile neutrino DM can be produced via the thermal freeze-in 
of s. In addition, h couples to N\v with suppressed coupling, so after the EW symmetry 
breaking there is a small mixing between and s. To check the effect of this mixing we 
consider the h —>■ u e Ni contribution to sterile neutrino production as well. The production 
by the singlet scalar has been considered in Ref. [21, 22], The thermal freeze-in production is 
caused by the Yukawa interaction of s and N\ or h and N\. Under the assumption m s m^, 
as the universe is expanding the temperature becomes low and s disappears first. The Higgs 
boson h, however, is still in thermal equilibrium and thermal freeze-in production by h is 
effective until T ~ m^. 

The dark matter yield can be calculated by solving the Boltzmann equations. In this 
scenario, the relevant Boltzmann equations for Y N] = /s are given as follows, 

dY Nl dYff dY™ 

(IT dT (IT ’ 1 ’ 

where Y^{Y^) is Yv, from s ( h ) decays respectively. 

In FIG. 1(a), FIG. 1(b) and FIG. 1(c), we show numerical results for the evolution of the 
sterile neutrino yield Yv, and the singlet scalar yield Y s for various thermal freeze-in mecha¬ 
nisms. Sterile neutrino DM is generated by the thermal freeze-in production of h in FIG. 1(a), 
the thermal freeze-in production of s in FIG. 1(b), and the non-thermal decay production of 
s in FIG. 1(c). 

Now, Ytv x is obtained from the following calculation. The Boltzmann equation for the 
sterile neutrino number density tin 1 involving s is written as, 

^ n Nl + 3 Hn Nl = ^ / d'K s d'K Nl d'K Nll (2'K) i 5 (i:) ( p s - p Nl - Pn^) 

spin 

{IVlLv.A p./.(1 - M (1 - /»'.) - \M\ 2 NiN ,„JnJn‘i (1 - /,)} • (12) 

We assume that the initial abundance of sterile neutrinos can be neglected and the singlet 
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(a) h -» v e Ni 


(b) s —> N\N\ 


(c) h —> s —> N\Ni 


FIG. 1. We describe the evolution of the yields Fv, and Y s as the temperature T decreases. The 
three thermal freeze-in production processes are shown in FIG. 1(a), 1(b) and 1(c). The sterile 
neutrino is generated by the thermal freeze-in mechanism of the Higgs boson (FIG. 1(a)), by the 
thermal freeze-in production of the singlet scalar (FIG. 1(b)), by non-thermal decay production of 
the singlet scalar (FIG. 1(c)). 


scalar enters into thermal equilibrium, such that, 

d 
dt 


n Nl + 3 Hn Nl = ^ / d^sdn^dn^^n^d^ (p s - p Nl - pjvq) \M\] feq 


f e 

\s—tNiN'i J s 


= 2 


spin 

I<i ( m a /T) 


r (s ->• N 1 N 1 ). 


(13) 


K 2 [m s /T)- 

The sterile neutrino yield Fv, = —p- can be obtained from the entropy density s = ^ h e ftT 3 , 
and satishes the following equation, 

1 1 Ki(m s /T) 


dYft 

dT 


45 


■K*G Ny /g^T*K 2 (m a /T) 
135a/5 m p 1 m 2 s K\ ( m s /T ) 


47^/2 h eSy /g^ T5 

Similarly, h also produces N\ as the following, 

dY& h = 1 1 IU (' m h /T) 

dT 


r (s -> N, Ah) Y s eq 
r ( s -> iViiVi). 


(14) 


4n 3 G N T3 K i ( m h/T) 
135a/5 rripi m 2 h I<i ( m h /T) 


r (h -> N x v e ) Y* q 


T (/r —> NiVe) , 


(15) 


87T 11 / 2 h eS T 5 
where m p 1 = 1.22 x 10 19 GeV is the planck mass, g e g and h e g are the effective number of 
degrees of freedom for energy and entropy and K n (x) is the modified Bessel function of the 
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second kind. The equilibrium yields Y^ q h are expressed as, 


Y eq = 4 5 g s m 2 s K 2 (m s /T) 
s 4vr 4 T 2 h eS ’ 
_ 45g fe m% /i 2 (rrih/T) 
h 4vr 4 T 2 /z eff 


(16) 

(17) 


The partial decay width of s into N\N\ is obtained as, 


T(s^N lNl ) = ^n- 


4M 


2 \ 3/2 


JVi 


m: 


and the partial decay width of h into N\v e is given as, 


r (h-> N { u e ) = 


Vel m h 

8tt 


M 


2 \ 3/2 
iVi 


mi 


n\m s 
167r 


y 2 eim h 

8tt 


(18) 


(19) 


In order to estimate the yield, we analytically integrate the relevant Boltzmann equations. 
The yield at the temperature of the universe today Y Nl (T 0 ) is given for s —>• iV 1 A r 1 as, 

y d * (t ) _ 135 ^ mpl [ To m 2 s K ^m s /T) )NN \ dT 

Y{To) - ~ 4X 11 / 2 h,„^iJ TR , ™ r (S ^ NM) dT 




135a/5 m p i m 2 s Ki ( m s /T ) 


47T 11 / 2 h eS , 




T (s —>■ NiNi) dT 


1.58 x 10 


14 l m N 1 


m a 


(S)J ’ 


( 20 ) 


where we assume h e s ~ (/ e ff ~ 100. The relic density of the sterile neutrino DM can be 
obtained as, 


9/^K 1 = 2.733 x 10 8 • l'o • 


fm DM \ 
V GeV / 


The DM relic density observed by Planck+WP [38] is estimated as, 


( 21 ) 


a DM h 2 = 0.1199 ± 0.0027. 


( 22 ) 


The sterile neutrino mass required to explain the observed DM relic density is thus m/v, « 
10 keV for (S) ~ m s ~ 1 TeV and my ~ 1 MeV for (S) ~ rn s « 100 TeV. FIG. 2 shows 
the relic density of sterile neutrinos as a function of m s for different values of my. 





















(a) s ->• ATiATi (5) = 10 TeV (b) s -> N X N ± (5) = 1000 TeV 

FIG. 2. The relic density of sterile neutrinos as a function of m s for different values of mjvj and 
(5) in the case of thermal freeze-in production via s. 


Similarly, we can integrate the Boltzmann equation of thermal freeze-in production via 
h, obtaining, 


r ATi 


(To) = 


135\/5 rripi f T ° m 2 h Ki (■ m h /T ) 


87r 11 / 2 T5 


r (/'i ->• iVn/ e ) dT 


135^/5 


rripi 


m\K i (m h /T) 


87T 11 / 2 h efr ^i 

1.04 x 10' • sin 2 6 ■ m 2 Nl , 


J'b 


r (/i ->• iVji/e) dT 


(23) 


with the relic density of N\ given as, 

Q^h 2 = 2.84 x lO“ 3 sin 2 0(^^) 3 . (24) 

Finally, sterile neutrino DM is also produced by the thermal background of active neutrinos 
via coherent scattering (Dodelson-Widrow mechanism). The dark matter relic density is 
found to be given as [39], 

9™K 2 = 5.47 x lO 7 sin 2 0(^V' 63 . (25) 

In the case of keV-MeV-scale sterile neutrino DM, the contribution to 9™ given in 
Eq.(25) is larger than that from the thermal freeze-in production via h given in Eq.(24). 
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In fact, there are additional contributions of the same order as those given in Eq.(24) that 
come from the decay of Z-bosons or W-bosons due to neutrino mixing, but we can safely 
ignore these contributions in the mass region under consideration. Altogether, the total 
sterile neutrino DM relic density is given as, 


jVi /i 2 


4.32 x 

VkeV. 
+5.47 x 10 7 sin 2 6 


3 / Te V \ 
\m s J 

rriN j _ \ 163 

.keV / 



(26) 


When the mass and the VEV of the singlet scalar are of order 1 TeV, the dominant mech¬ 
anism for production of keV-MeV-scale sterile neutrinos is via the thermal freeze-in of the 
singlet scalar. 


B. The singlet scalar is out of thermal equilibrium 


In this case A <C 10 -6 and both s and N\ never enter into thermal equilibrium in the early 
universe. Sterile neutrino DM is generated via the thermal freeze-in of h. The singlet scalar 
s is also generated by thermal freeze-in production, and then proceeds to decay efficiently 
into N\. In this section, we assume rn s <C mjv 2 , 3 , so that s cannot decay into N 2 and N 3 . 

To calculate the yield of sterile neutrinos, we have to solve the Boltzmann equations given 
by the following two interrelated equations, 

dY s dY A dY° s 

~dT ~ ~dT~ dT ’ 

dY Nl dY° s dY° h 

dT dT + dT ’ 


(27) 


d Y 


Ds 


dT 


1 dY°‘ 

2 dT ' 


(28) 


Standard model particles in thermal equilibrium can annihilate into a singlet scalar. 
For simplicity, we concentrate on thermal Higgs annihilation as the dominant production 
mechanism of singlet scalars and ignore the other standard model effects. The Boltzmann 
equation for the annihilation process can be expressed as follows, 


dY' 


135\/5 


m p i 


A 2 m 6 


dT 64 tt 17 /2 h effv ^i T 3 


K x (2 m s /T). 


(29) 


We integrate Eq.(27) to estimate the yield of sterile neutrinos Y^ x . There is no initial yield 
(Y s (T re ) = 0) and no final yield (Y s (T 0 ) = 0), and therefore, the following equation can be 
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(a) h —> s —» NiNi A = 10 7 (b) h —> s —> IVi-M = 10 keV 

FIG. 3. The relic density of sterile neutrinos as a function of m s for (a) different sterile neutrino 
masses mjvj and (b) different values of the Higgs portal coupling A. 


obtained, 


r-To dY Ds 

71 dT 


dT = 


pT ° dy/ 

"dT" 


dT. 


(30) 


The yield of sterile neutrinos at today’s temperature can be obtained using Eq.(29) and 
Eq.(30), 


Y& (To) = -2 


rT 0 jyDs 

s -dT = 2 


dT 


d s -dT 


Jt re 

135\/5 m p i 


JTre 

rT ° A 2 m 


dT 


327T 17 / 2 /r effv ^i T 3 

Therefore, the yield at today’s temperature Tv, (T 0 ) is given as, 

135\/5 m p i 


K\ (2m s /T) dT. 


(31) 


Yn : : (To) = 


r ‘ T ° A 2 m< 


327T 17 / 2 deffi/deff J Tre T 3 
135\/5 m p i 7° A 2 m s 


Ad (2m s /T) dT 


327T 17 / 2 


T 3 


ATi (2m s /T) dT 


A 2 

« 1.07 x 10 13 ( — 

m s 


(32) 


The sterile neutrino DM relic density resulting from the non-thermal decay mechanism is 
obtained as, 
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FIG. 3 shows the relic density of sterile neutrinos as a function of m s for different values of 
the Higgs portal coupling A and the sterile neutrino mass m n 1 . 

In this scenario, the Dodelson-Widrow mechanism can also produce sterile neutrino DM. 
Therefore, the total relic density is obtained as, 

,2 


O(p)ly) 


+5.47 x 10 7 sin 2 #(^) 
Vkev / 


1.63 


(34) 


The relic density formula depends on the Higgs portal coupling A. The coupling A is bounded 
as A < 10 -6 so that s does not come into thermal equilibrium. Therefore, s can not be 
heavier than the TeV scale in order to produce keV-MeV sterile neutrino dark matter via 
the non-thermal decay production mechanism. 


IV. X-RAY CONSTRAINTS AND LIFETIME BOUNDS ON STERILE NEU¬ 
TRINO DARK MATTER 

In this section, we will review the X-ray bounds and the lifetime bounds on sterile neutrino 
DM. Sterile neutrinos can decay into standard model particles through active-sterile neutrino 
mixings. In the keV-MeV mass range, sterile neutrinos decay mainly into the three active 
neutrinos [40-42]. For the three-neutrino decay channel, the decay lifetime is expressed as, 

iq /keVA 5 1 , , 

t 3u ~ 2.88 x 10 19 sec - —35 

\m Nl J sin' 0 

Their lifetime must be longer than the age of the universe (10 17 sec) if sterile neutrinos are 
to constitute dark matter, which constrains the mixing angle and the sterile neutrino mass 
as follows, 

sin 2 20 < 2.88 x 10 2 f^fk) 5 . (36) 

If sterile neutrinos constitute dark matter, their radiative decay (TV] - 471 /) would lead 
to a cosmic X-ray background. We have not seen such an X-ray excess except for the recent 
observation of a 3.5 keV signal [43, 44] in galactic clusters. This puts an upper limit on the 
neutrino mixing angle for a given sterile neutrino mass. From the diffuse X-ray background 
observations XMM-Newton [45, 46] and HEAO-1 [47], the authors of Ref. [19, 48] obtain the 
simple empirical formula, 

sin 2 20 < 1.15 x . (37) 

V keV / VO-26 y V ; 
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The XMM-Newton observations of the Virgo and Coma galaxy clusters present the more 
stringent constraints [19, 49], 



(38) 


More precise X-ray constraints have been reported in Ref. [19]. Note that these bounds are 
given for sterile neutrino DM which explains the current dark matter density. If sterile 
neutrino DM only constitutes part of the total dark matter, the X-rays bounds become 
weaker. 

V. THE FREE STREAMING HORIZON AND LYMAN-a CONSTRAINTS 

Recent observations such as the WMAP and Planck missions have proven that the ACDM 
model, which contains cold dark matter, is an extremely successful cosmological model [38] . 
However the ACDM can not solve the small-scale crises [50], including the missing satellite 
problem and the cuspy halo problem. Warm dark matter (WDM), which has an adequate 
free streaming horizon and suppresses the structure of dwarf galaxies size, may solve the 
problem. The upper bound on the free streaming scale of WDM is obtained from the ob¬ 
served Lyman-a forest, which refers to the absorption lines of intergalactic neutral hydrogen 
in the spectra of distant quasars and galaxies. 

The free-streaming horizon corresponds to the average distance travelled by DM particles 
and is a good measure to classify CDM, WDM and HDM. The free streaming horizon is 
given as, 



(39) 


where t m is the DM production time, to is the current time, {v (t)) is the average thermal 
velocity of the DM particles, and a(t) is the scale factor. In this paper we assume that 
the free-streaming scale of CDM, WDM and HDM satisfy A fs < 0.01 Mpc, 0.01 Mpc < 
Aps <0.1 Mpc and 0.1 Mpc < Aps, respectively. This is not an accurate definition, but 
gives a useful criteria to classify the thermal property of DM. Note that HDM is excluded 
by observations of the Lyman-a forest. 

In order to determine the free-streaming horizon, we now consider the average thermal 
velocity (v (t)) of the sterile neutrino DM N { . We define t nr as the time when Ni becomes 
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non-relativistic, which we take to be when the equality ( p(t nr )) = rri /v, is satisfied. The 
approximate average thermal velocity (v ( t )) is then given as follows, 


{v (t)> ~ 


1 t t nr . 


(pft)) 

m N , 


t ^ t"nr • 


(40) 


The non-relativistic thermal velocity is expressed in terms of the average thermal momen¬ 
tum, which can be extracted from the distribution function / (p) and depends on the DM 
production mechanism. In this section we will consider the average thermal momentum and 
the free streaming horizon when production is via thermal freeze-in of the singlet scalar, 
via the Dodelson-Widrow mechanism and via the non-thermal singlet scalar. Finally, we 
determine the Lyman-a constraints and the allowed parameter region for each production 
mechanisms. 


A. Production via thermal freeze-in of the singlet scalar 

For production via the thermal freeze-in of the singlet scalar boson, the momentum 
distribution of sterile neutrino DM [51, 52] is given by, 


where 


/ M = aw ( p/t) 


9u o) = y 


n= 1 


rr 


(41) 


(42) 


The normalization factor (3 is determined by the Yukawa coupling n and the singlet scalar 
mass m s , with j3 oc . The average thermal momentum ( p{t )) can be calculated as, 

5? dp E“! 


(p (*)> = 


,5/2 




2.4527 T. 


(43) 


.5/2 


This average thermal momentum {pit)) leads to the average thermal velocity (v (t)), 


(v {t)) ~ 


2.45T _ dl(^nr) 

a(t) 


m N . 


t ^ tnr-) 


t>t r 


The time when DM particles become non-relativistic is th( 2 « 2.45 sec. 


(44) 


The free 
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(a) s -> TViTVi <5) = 1 TeV 


(b) s -> N\Ni ( S ) = 100 TeV 


FIG. 4. X-ray bounds, the free streaming horizon, HDM, WDM and CDM regions and other 
constraints for sterile neutrino production via thermal freeze-in of the singlet scalar for (a) (S) = 
1 TeV and for (b) (S) = 100 TeV. 


streaming horizon is calculated as, 

(v («)> 


X F s — 


a (t) 
dt 
a (t) 


dt 


+ 


{v (*)) 

a (■ t ) 


rto 


dt T 


dt 


5y/t e qt. 
® (teq) 

\Jteqtnr 

® (^eq) 


nr y/teqt 


® ideq) 

5 + In ( — 

trt.r 


nr In f ^ ^y/^eqt'. 


a (■ t ) 

^^teqtr 
a (teq) a (teg) 


1/2 


(45) 


To obtain the last line, we neglect the third and the last terms of the third line. 

In this production mechanism, the DM is produced at high temperatures, T > 1 TeV, 
and entropy dilution affects the free streaming horizon. The effect of entropy dilution can 
be estimated by the factor which is given by, 


g eS (high T) _ 109.5 
g e s (current T 0 ) 3.36 


(46) 


Now we assume that both s and N\ contribute to the effective number of degrees of freedom 
and ignore the tiny effect of the other heavy right-handed neutrinos N 2 and N ?J . Taking 
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entropy dilution into account and using the conversion factor c = 10 14 (Mpc/sec), the final 
expression is given as, 



(47) 


The Lyman-cx bound on is given by, 


m/v, > 1.57 keV. 


(48) 


The range of sterile neutrino mass corrsponding to WDM is obtained as, 


1.57 keV < m Nl < 20.5 keV. 


(49) 


In FIG. 4, we show the X-ray bounds and the HDM, WDM and CDM regions for sterile 
neutrino production via the thermal freeze-in of the singlet scalar. In this figure, we assume 
that m s is larger than m h = 125 GeV but smaller than (a) ( S) = 1 TeV and (b) ( S) = 
100 TeV. We also show the parameter region where more than 1 % of the DM is produced 
by the DW mechanism. When ( S ) is 1 TeV (100 TeV), the sterile neutrino DM is warm 
(cold). However, the scenario with ( S) > 100 TeV suffers from the X-ray constraints. 

B. Production via non-thermal decay of the singlet scalar 

If the Higgs portal coupling is small and the singlet scalar is out of thermal equilibrium, 
it decays into the sterile neutrino. The free streaming horizon was considered in [23, 24], 
The momentum distribution of the sterile neutrino DM is given as [53-56], 



(50) 


where f3 is a normalization factor and the DM temperature is T^m — • ^ ie avera g e 

thermal momentum is given as, 


From Eq.(51), the average thermal velocity is expressed as, 



(51) 


1 




(52) 


y/7rm s a(ti n ) _ a(t 

nr') 

4 m^ 1 a(t) a(t) 
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FIG. 5. X-ray bounds, HDM, WDM and CDM regions and other constraints on the sterile neutrino 
mass mjVi and mixing angle 6 in the case of production via the non-thermal decay of the singlet 
scalar for (a) rn s = 1 TeV and for (b) rn s = 100 TeV. 

Now, we assume that the production time is ti n = tf e + r, where tf e is the freeze-in time of 
s and is given as tf e ~ sec > w h ere Tf e ~ m s GeV is the freeze-in temperature, and 

the lifetime of s is r = h/T (s —> N\ N \). The time at which the sterile neutrinos become 
non-relativistic is given by t nr = f-U n sec and the time of matter-radiation equality is 
t eq = 1.9 x 10 11 sec. We estimate the free-streaming horizon of the DM sterile neutrinos 
using the formula, 


Xfs — 


C\/teqtr 
a (t. 


eqj 


5 + In ( — 

tn.r 


e /r 


For m s = 1 TeV and K\ = 10~ 8 , the Lyman-a bound on m/v, is obtained as, 

m at, > 4.36 keV. 

The WDM sterile neutrino mass can be constrained as 

4.36 keV < m/v, < 64.3 keV. 

For m s = 100 TeV and K\ = 10~ 8 , the Lyman-a bound is given as, 

m/v, > 64.2 keV. 


(53) 


(54) 


(55) 


(56) 
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The WDM sterile neutrino mass range is obtained as 


64.2 keV < m Nl < 840 keV. 


(57) 


Therefore, in this scenario, the singlet scalar can not be heavier than the TeV scale 4 This 
constraint is tighter than that in the DW mechanism m^ 1 > 10 keV [57, 58]. This is because 
the DM sterile neutrino is produced by the decay of non-thermal heavy particles. When the 
lifetime of the singlet scalar is not small, the Lyman-ct constraint becomes tight. 

In FIG. 5, we show the X-ray bounds and the HDM, WDM and CDM regions in the 
rrii\T l -6 plane for sterile neutrino DM production via non-thermal decay of the singlet scalar. 
In FIG. 5(a) we take m s = 1 TeV and plot the sterile neutrino constraints for the three 
different values of Higgs portal coupling A = 10“'' 4 ,10“' 7 ,10“ 8 . The sterile neutrino DM is 
not constrained by the Lyman-ct bounds. FIG. 5(b) shows the sterile neutrino constraints for 
the portal couplings A = 10“ 6 ' 4 ,10“ 6 ' 7 ,10“ 7 and singlet scalar mass m s = 100 TeV. When 
the mass of the singlet scalar is large, the produced sterile neutrinos are warmer and the 
scenario is constrained by the Lyman-a bounds. 

C. Sterile neutrino DM production mechanisms in the z/MSM 

In the i/MSM, sterile neutrino dark matter can be generated by the thermal freeze-in 
mechanism - through interactions with the Higgs boson - or the Dodelson-Widrow mech¬ 
anism. The free-streaming horizon of sterile neutrino DM produced by thermal freeze-in 
via the Higgs boson is as small as that of sterile neutrino DM that is produced by thermal 
freeze-in via s. However, this production mechanism is in conflict with lifetime bounds and 
X-ray bounds (see FIG. 6(a)). 

In the Dodelson-Widrow mechanism, sterile neutrinos are produced from the thermal 
background of active neutrinos via coherent scattering. Therefore, the momentum distribu¬ 
tion is thermal and of the Fermi-Dirac type [14], i.e, 



(58) 


where p denotes the comoving momentum of N\ and (3 oc 6 2 M\. For the thermalized sterile 


4 Ref. [27] presents more detailed calculations in this scenario, solving numerically the system of Boltzmann 
equations. The singlet scalar mass could be more tightly restricted. 
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the parameters for which the sterile neutrino DM density QdmP 2 = 0.1199 and the red dashed 
lines the parameters for which the sterile neutrino DM density satisfies CI^mP 2 = 0.01199. 


neutrinos, the average thermal momentum (p (t)) is given as, 

7vr 4 T 


f°° p 

/ f-t-W ^0 ^e p / T -\-l 

(PW) = ~roo~, -7^“ 

Jo d P~e 


180C(3) 


3.1513 T. 


(59) 


io ^e eP /T + 1 

The average thermal momentum of the thermally produced sterile neutrinos thus satisfies 
the relation (p) /3.15 T ps 1 , but the DW mechanism produces a colder distribution of sterile 
neutrinos, with (p) /3.15 T « 0.9. The free-streaming horizon of the sterile neutrino DM is 
given in terms of the average thermal momentum as [41], 


\ FS ps 0.84 Mpc 


keV 

m Nl 


(P) 


3.15 T 


(60) 


In the case of the DW mechanism, observations of the Lyman-a forest lead to the severe 
constraint rn^ l > 10 keV [57, 58]. This limit is in conflict with X-ray bounds (see FIG. 6 (b)), 
meaning that the DW mechanism scenario is excluded. 

If the lepton asymmetry is relatively large in the early universe, the thermal production 
of sterile neutrinos can be enhanced by the MSW effect (Shi-Fuller mechanism). The Shi- 
Fuller mechanism leads to a colder thermal distribution with (p) /3.15 T ps 0.6. Therefore, 
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the Shi-Fuller production mechanism can evade the Lyman-ct bounds and X-ray constraints, 
although a relatively large lepton asymmetry is needed. In conclusion, the thermal back¬ 
ground production of sterile neutrino DM, through mechanisms such as the DW mechanism, 
is severely constrained due to the large free-streaming scale and X-ray bounds. It should be 
emphasized that the z/MSM fails to explain the dark matter sector, because the sterile neu¬ 
trino dark matter cannot be generated by neither the DW mechanism nor thermal freeze-in 
production via the Higgs boson. This is why we have to extend the z/MSM to explore new 
DM production scenarios. 


VI. THERMAL FREEZE-IN LEPTOGENESIS VIA THE SINGLET SCALAR 


In this section, we discuss leptogenesis scenarios that rely on the thermal freeze-in produc¬ 
tion mechanism. In this scenario, we assume that the heavy Majorana neutrinos are gener¬ 
ated by either thermal freeze-in via the singlet scalar or non-thermal decay of the frozen-in 
singlet scalar. The produced Majorana neutrinos generate a lepton asymmetry, which is 
transferred into a baryon asymmetry via non-perturbative electroweak effects (Sphalerons). 
We assume that N { constitutes sterile neutrino DM and does not affect the leptogenesis 
scenarios, while N 2 and N 3 satisfy Tew < M 2 < M 3 . 

In line with the original motivation for the z/MSM to explain BSM phenomena by intro¬ 
ducing TeV-scale particles, we set the Majorana masses to around the TeV scale. However, 
TeV-scale leptogenesis is in conflict with the Davidson-Ibarra bound [59], which constrains 
the Majorana mass as M 2 > 10 9 GeV. It is possible to evade this lower bound when the 
mass difference between N 2 and IV 3 is of the order of their decay width. In this case, resonant 
leptogenesis can occur [60, 61]. The resonant CP asymmetry is obtained as, 

V(N i ^£ a H)-V{N l ^T a H*) 


e* = 


« / / - 
~r. -r-r \ ~ ^ ^ ^ 


T (Ni -a l a H) + T (Ni -a £ a H*) 

The e'-type CP asymmetry is obtained from the vertex contribution, 


(61) 


£; = 


Im {y ] y) 
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Ni 


>33 


/ 


m 


Ni 


m 


Ni 


(y*y)ii(y*y) 

where T ^ is the tree-level decay width and f{x) is the loop function, which are given as, 

T + x 


(62) 


(y ] y ),, , x 

r V = o rn N f (x) = 


87T 


(1 + x ) In 


x 


(63) 
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In the degenerate heavy Majorna neutrino mass limit (mjv 2 ~ i^n 3 ), the CP asymmetry e' 2 
cancels with e' 3 and no net CP asymmetry can be obtained from the vertex contribution. 
The £-type CP asymmetry arises from the self-energy contribution, 


£i = 




{y ] y) ii {y ] y) jj \m N . 


Ni 


Vi - m N,) m Ni m Nj 


+ m N? ri 


(64) 


In the limit (mjv 2 
In the limit m^ 2 


s mjv 3 ), £% dominates over Furthermore, e 2 and £3 have the same sign, 
mjv 3 , the total CP asymmetry involving the N 2 contribution is given by, 


£2 - 


Im {y ] y) 


2 

23 


Wv)mWy )33 


N 3 


( m N 2 ~ m N 3 ) m N 2‘m N3 


(mL - m 2 r 


Y + 


°N 2 "°n 3 j P 111 X 2 ‘ x 3 

In order to generate an 0(1) lepton asymmetry, it is necessary to satisfy the following two 
conditions [60], 


2 p2 


(65) 


?7Tjv 3 — 77lj\r 2 


2 1 ^3,2) 


1 . 


(j/ty) 22 (2/ + 2/) 33 ~ ^ 

In general, the mass difference is larger than the tree-level decay width (Amjv 32 > T Ns2 ), 

such that, 

Im (z/ + y) 23 t m N 2 - 171 N 3 ) m N 2 m N3 


C 2 - 


(m ^ 2 - m ^ 3 ) 2 + m^ 2 r ^ 3 
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When we assume 


(2/ t 2/) 22 (2/ t 2/)33 Am ^32 
Im (y + y) 2 3 _ 10 _ 3 


) 2 

_ / 23 

M 22 M 3 


mjv 2iS = 1 TeV and ArriN 32 


(67) 

1 MeV, the CP 


asymmetry factor becomes £2 ~ 10 ~ 9 . 

In thermal leptogensis, the right-handed neutrino yeild Yjv 2 and the lepton asymmetry 
Yal satisfy the following two Boltzmann equations, 

dY N2 


= (D 2 + S) ( Y n ■„ - Yg) . 
rjT = -c 2 D 2 (Y n , - Yg) + W, D Y AL . 


dT 

dV AL 


( 68 ) 

(69) 


The scattering term S describes A L = 1 scattering effects, but we neglect this contribution 
for simplicity. The decay and washout terms are expressed as, 

r iv 2 1 Ki (■ m^/T) 


D ^ = im¥i ^ 2 wy w ’ d ^ = \ d ^ t) YJA w^ 


(70) 
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The equilibrium yields of Majorana neutrinos and leptons are given by, 

\/-eq (rp\ _ 45?7 4 K 2 ( m N jT) ^ eq _ 15 ^ 

Y nA 1 )- 2^4 T 2 heS ’ Y i - 4n 2 heS ^ 

The analytical solution for the lepton asymmetry Yal is given by the following formula [62- 
64], 

Y al (T) = Y al (The) e S ^ E dT ' w ^ T ') _ f dT'e 2 ^e# dT " w ^ T ") (72) 

Jt re 


— -62 


dT , dYjy 2 JT dT n WiD (T ») 

(IT 


(73) 


If we assume that there is no preexisting lepton asymmetry, Yal (Tre) = 0, and neglect the 
washout term, Eq (73) reduces to the form, 


f T o f iy 

Yal (To) ~ -e 2 / dl* = e 2 Y Na (Tre) 
Jt re 


(74) 


In general, the heavy Majorana neutrinos come into thermal equilibrium. The initial yield is 
given by the thermal equilibrium yield Yv 2 (Tre) = Y^ ~ 0.004 and the lepton asymmetry 
is approximately given by Yal (To) — 0.004 e 2 . If we include washout effects, the final lepton 
asymmetry can be obtained as [64], 


Yal(To)^~€ 2 ( Fn2 


Y%. 


(75) 
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So far we have been describing the thermal leptogenesis scenario without the singlet 
scalar. Now we discuss how this scenario will be modified by the singlet scalar, both when 
it is in and out of thermal equilibrium. 


A. Leptogenesis via the singlet scalar in thermal equilibrium 


In this subsection, we discuss a new leptogenesis scenario where the singlet scalar is in 
thermal equilibrium and the Majorana neutrinos are generated by thermal freeze-in via the 
singlet scalar. In this scenario, the Higgs portal coupling has to be relatively large, A > 10~ 6 , 
and the Yukawa coupling needs to be small, k 2 < 10~ 6 . The relevant Boltzmann equations 
are given by the following formulas, 

dY N2 


dT 

dY AL 

dT 


= D 2 (Y n , - Yg) - 2 D.Y?, 

= -e 2 B 2 (IV - iV) + W, d Y al . 
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The relevant terms are given by, 
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D s Y s eq = - 
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T (s -> 7V 2 7V 2 ) 
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We can write down the analytical solution for the lepton asymmetry Yal as, 

Y al (T) = -e 2 f T + 2A s Y/ 9 ) e-lr' dT ” w " } ^dr. 

JTre V / 


(80) 


(81) 


If we assume that the Majorana neutrinos and the initial lepton asymmetry are zero, Yv, 2 3 = 
Y A l (Tre) = 0 , and neglect the washout term in Eq.(81), then we find, 


Yal (To) ~ 6 2 Ytv 2 (T rb ) - e 2 f ° 2A S Y / 9 dT 

•'Tre 

rTo 


- -^2 


2£> S Y / 9 (IT. 


(82) 

(83) 


Finally, we analytically integrate Eq.(83) from To = 0 to Tre = oo and obtain the lepton 
asymmetry as, 
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(84) 


The (A + L)-violating interactions of sphalerons come into thermal equilibrium at tempera¬ 
tures above the electroweak phase transition T > T c & 200 GeV, and the lepton asymmetry 
can be converted into a baryon asymmetry as follows [65-67], 


OO OO 

Tab (T) = -Yab-l (T) = --Yal (T). 


(85) 
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(a) s —> N 2 N 2 e 2 = 10 9 (b) s —> N 2 N 2 k 2 = 10 7 

FIG. 7. These figures show the dependence of the baryon asymmetry yield Tab on the CP asym¬ 
metry factor 62 and Yukawa coupling k 2 ■ The Yukawa coupling should satisfy n 2 < 10~ 6 , otherwise 
the sterile neutrino is in thermal equilibrium. 


Therefore, the final baryon asymmetry is given by, 
Yab (To) « 0.87 x 10~ 12 ( 62 
« 0.87 x 10 2 


10-9/ V10- 7 
62 


k 2 ^ ^ / TeV\ 

\m s J 


10 
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) \TeVJ V m 8 J V (S') 


From BBN results, the baryon asymmetry is determined to be, 


-li 


y ab = (8-!0 ± 0.85) x 10 


From CMB measurements, the baryon asymmetry is determined to be, 


( 86 ) 


(87) 


Y2b B = (8.79 ± 0.44) x 10 


-11 


( 88 ) 


FIG. 7 shows the dependence of the baryon asymmetry on the CP asymmetry factor 62 and 
the Yukawa coupling k 2 . The constraint on the Yukawa coupling k 2 < 10 ~ 6 and Eq.( 86 ) 
lead to the following constraint on m s , 


m s \ 
TeV/ 


< 1.09 


62 

IQ " 9 


(89) 
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Therefore, the mass of the singlet scalar cannot be larger than 1 TeV if we are to produce the 
observed amount of baryon asymmetry. In section III A, we discussed the thermal freeze-in 
production of keV-MeV sterile neutrino DM and concluded that the singlet scalar should not 
be heavier than the TeV scale. The mass of the singlet scalar required to achieve leptogenesis 
is more severely restricted than the mass required for successful dark matter scenarios 5 . 


B. Leptogenesis via the singlet scalar out of thermal equilibrium 


When the Higgs portal coupling is small, A -C 10“ 6 , and s, N 2 and A 3 do not exist 
in the early Universe, they do not come into thermal equilibrium. The singlet scalar is 
then produced by the thermal freeze-iu mechanism and decays efficiently into the Majorana 
neutrinos N 2} 3 which generate a net lepton asymmetry. In this scenario, to determine the 

yields we have to solve the following Boltzmann equations, 

dY, dY Q A 


dT 

dY]sr 2 

dT 

dY AL 

dT 


dT 


+ d s y s , 


= D 2 (Y n . 2 - y£) - 2 D s Y a , 

= 6 2 D 2 (Y N2 - Yg) + W id Y al . 


The relevant term in the annihilation process can be expressed as, 

dY A 135a/5 7T7.pl A 2 777 . s 
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dT 64tt 17 /2 heS T 3 

The singlet scalar decays into the Majorana neutrinos, and the decay term is given by, 

1 1 Ad (■ m s /T) 


(90) 

(91) 

(92) 

(93) 


D,Y, = 


45 


47r 3 Gjv T 3 K 2 ( m s /T) 
3a/ 5 777 p i 1 K \ (■ m s / T ) 


T (s -)• N 2 N 2 ) Y s 


r (s —>■ n 2 n 2 ) y s . 


2tt 3 / 2 ^/g^g T 3 Ii 2 (m s /T) 

We then integrate Eq.(90) to estimate the yield of Majorana neutrinos Vv 2 , 
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Taking the initial yield Y s (T RE ) and the final yield Y s (T 0 ) to be zero, the following equation 
is obtained, 


r-To 


BY AT = - 


rTo dYd 


'Tre 
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dT 
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(96) 


5 From Eq.(21) and Eq.(86) we see that in order to satisfy observational constraints on the DM density and 
baryon asymmetry in this scenario, the mass of the sterile neutrino must be around 100 MeV. 
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(a) h —> s —> N 2 N 2 e 2 = 10 9 (b) h —> s —» N 2 N 2 n 2 = 10 ' 

FIG. 8. These figures show the dependence of the baryon asymmetry yield Iab, on the CP 
asymmetry factor e 2 and the Higgs portal coupling A. The Higgs portal coupling has to be small, 
A < 10 -6 , in order to prevent the singlet scalar from entering into thermal equilibrium. 


The lepton asymmetry Yal at the temperature of the universe today is obtained as, 
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The final baryon asymmetry is expressed as, 
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(98) 


In this scenario, a TeV-scale singlet scalar can generate the observed amount of baryon 
asymmetry. FIG. 8 shows the dependence of baryon asymmetry on the CP asymmetry factor 
e 2 and the Higgs portal coupling A. The constraint on the Higgs portal coupling A < 10~ 6 
and Eq.(98) lead to the following constraint on m s , 

^< 0 . 13 ^ 62 
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In this leptogenesis scenario, the singlet scalar cannot be heavier than 1 TeV in order to 
produce the observed baryon asymmetry. In section III B, we considered the production 
of keV-MeV sterile neutrino DM from the non-thermal decay of the scalar singlet, and we 
showed that the singlet scalar needs to be lighter than the TeV scale. This leptogenesis sce¬ 
nario constrains the mass of the singlet scalar more severely than the scenario that explains 
the sterile neutrino DM 6 


VII. CONCLUSION 

In this paper we have considered an extended z/MSM with one additional singlet scalar. 
The existence of the singlet scalar and the associated freeze-in production of sterile neutrinos 
and Majorana neutrinos make it possible to alleviate the usual tensions found when trying 
to construct dark matter and leptogenesis scenarios within the framework of the z/MSM. 
We have studied two scenarios, making different assumptions about the thermal properties 
of the singlet scalar in each scenario. If the Higgs portal coupling is relatively large, so 
that the singlet scalar enters into thermal equilibrium, then the sterile neutrino and the 
heavy Majorana neutrinos are produced directly via the thermal freeze-in mechanism. If, 
on the other hand, the Higgs portal coupling is much smaller and the singlet scalar is out of 
thermal equilibrium, first the singlet scalar is produced by the thermal freeze-in mechanism. 
Then, the sterile neutrino and heavy Majorana neutrinos are produced via the non-thermal 
decay of the singlet scalar. In these scenarios, the sterile neutrino DM can evade Lyman-a 
bounds and X-ray constraints. We found the latter scenario to be more tightly constrained 
by Lyman-a bounds, with the singlet scalar needing to be lighter than 100 TeV to produce 
keV-MeV sterile neutrino DM. Thermal freeze-in leptogenesis scenarios severely restrict the 
singlet scalar mass to be less than 1 TeV in order to generate the observed baryon asymmetry. 
In summary, in the extended z/MSM with a singlet scalar, a TeV-scale singlet scalar mass is 
needed in order to generate the observed abundance of dark matter and baryon asymmetry 
via the thermal freeze-in production mechanism. 

6 Note that one singlet scalar cannot be used to explain both the dark matter sector and leptogenesis if we 
consider thermal freeze-in production of the sterile neutrinos. This is because the production of heavy 
Majorana neutrinos via non-thermal decay of the singlet scalar dominates. 
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